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Abstract. This paper was written to postgraduate students of Physics. It is 
shown didactically how to calculate rigorously the wrinkling wavelengths 
observed when rigid thin films are deposited on soft elastic substrates. Our 
predictions are compared with DLC films deposited on PDMS substrates. 

Key words: wrinkling wavelengths; thin films. 


(I) Introduction. 

In basic courses of Physics and Engineering we learned that it is not 
enough only know the forces applied on a structure. [1] It is also necessary to 
know the tensions and the deformations that they produce. All real 
substances are deformed under the influence of a forces and tensions. The 
part of Physics that studies these changes is usually known as "Theory of 
Elasticity. " [2,3] As this theme has been extensively studied in many books 
and papers we will make a brief review only of basic concepts adopted in 
the theory of elasticity. The mechanics of solid bodies, regarded as 
continuous media, forms the content of the Theory of Elasticity. [1,2,3] 

So, let us consider a bar of length i 0 with cross section area A 
submitted to a force F applied perpendicularly over the area A. We define 
stress o = F/A (force / area). If, under the action of force, the bar undergoes 
a longitudinal strain s c (strain) (compression or distension) given by 
= IAf| H 0 = |f - f 0 |/f G we define modulus of elasticity or Young 
modulus by 


E = o/(Af/f 0 ) = o/s£ (1.1). 

If the bar undergoes a longitudinal strain £ ( = Af/f o and a transverse 
strain s t = Ae/e 0 the relationship between the two deformations is defined as 

"Poisson's coefficient" (v): [4] 

v = - (Ae/e 0 ) / (Af/f 0 ) = - s t /s c (1.2). 
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If a body is deformed from an angle cp due to a shear stress o we 
define shear strain by = Ax/f G . The shear modulus of the material due 
to the strain g^ is defined by G = p: 

G = p = o/(Ax/f 0 ) = o/gq, (1.3). 

The shear modulus is also known as Stiffness modulus or torsion 
modulus. 

A body with volume V G submitted to a hydrostatic pressure P suffers 

a volumetric strain g v = AV/V 0 = (V-V 0 )/V 0 . The volumetric modulus 

K or bulk modulus (B) is defined by 

K = B = - P/(AV/Vo) = - P/g v (1.4). 

The inverse of the bulk modulus (K = B) is called compressibility k : 

k = 1/B = 1/K (1.5). 

One can show [4] the following Lame relations between B, G, E and v: 

K = B = E/3(l-2v) (1.6). 

and 

p = G = E/2(l+v) (1.7). 

In reference [4] are seen Lame relations between all elasticity coefficients. 


(1) Element of Volume Submitted to an External Force F. 

Under the action of applied forces, the solid bodies exhibit 
deformation to some extent, i.e. they change in shape and volume. The 
deformation of a body is described mathematically as follows. The position 
of any point P in the body is defined by its position vector r = (xi, x 2 , x 3 ) in 
some coordinate system. When the body is deformed, every point is in 
general displayed. The position of P after the deformation will be r'(with 
coordinates x3). The displacement of this point P due to the deformation is 
given by the vector u = r'- r, called displacement vector which is a given 
function of x i? that is, u = u(xi). This means that the coordinates of the 
displaced point P are functions of the coordinates Xi. [3] 

When a body is deformed the distances between its points change. 
Let us consider two points very close together. If before the deformation 
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we have distances dxj these distances after the deformation would be given 
by dx'i = dxj + du;. The original distance between the points df = Vdx ; 2 after 
the deformation would be given by df = Vdx', 2 . Since dx/ = dx; + dui and 
that dui = (5ui/5x k )dx k we get 

df' 2 = df 2 + 2u ik dx t dx k (1.1), 

where u ik is a symmetric tensor named strain tensor given by 

u ik = (l/2)[ (<3uj/ox k ) + (ou k /oXj) + (5uc/5xi)(0u c /5x k )] (1.2). 

Let us consider a cubic element of volume AV = Ax Ay Az of a 
body. When submitted to an external force it creates stresses Gy along the 
cube surfaces seen in Figure 1 generating tensor strains u* given by (1.2). 



Figure 1 . The strain tensor on the cube surfaces. The resulting moment of the shear 
forces must vanish'^ 

When u = u(x,y,z) are very small deformations products like 
(5ui/5x k )(5u k /5xi) can be neglected and (1.2) can be written as 

u ik = (1 /2) {du-Jdx k ) + (ou k /ox,)} (1.3). 

Any deformation can be represented as the sum of a pure shear and a 
hydrostatic compression. To do so we need only to use the identity [3] 

u ik = [u ik - (1/3)8* u e( ] + (1/3) 8 ik u C£ (1.4). 


The first term on the right is evidently a pure shear, since the sum of its 
diagonal terms is zero (Sn = 0). The second is a hydrostatic compression. 
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(2) Force F on the Body. 

Thus, for any portion of the body, each of the three force 
components |Fi dV of the resultant of all the internal stresses can be 
transformed into an integral over the surface. As is known from vector 
analysis, the integral of a scalar over an arbitrary volume can be 
transformed into an integral over the surface if the scalar is a divergent of a 
vector. In our present case we have integral of a vector, and not of a scalar. 
Hence, the vector F; must be the divergent of a tensor of rank 2, that is [3, ] 

Fi = do ik /dx k (2.1). 

So, the force on any volume of the body can be written as an integral over a 
closed surface bounding that volume: [3] 

Jv Fi dV = I v (do ik /dx k )dV = j s o ik da k (2.2), 

where o ik da k is the i th component of the force on the surface element da. 


(3) Free Energy of the Thermodynamic Deformation. 

Assuming valid Hooke's law, a general expression according to 
Thermodynamics of the free energy F per unit of volume of a deformed 
isotropic body is obtained summing two independent squared scalars of 
two components: one due to pure shear and another due to pure 
hydrostatic compression shown in (1.4): 

F = n [Uik - (l/3)5 ik u„ ] 2 + (1/2) K u tt 2 (3.1). 

It can be also shown that, 3 

o ik = dF/du ik (3.2), 

from which one can determine the stress tensor o ik : 

G ik = 2|i [u ik - (1/3)8* u a ] + K u u 5 ik (3.3). 

This expression determines the stress tensor a* in terms of the strain 
tensor u ik for an isotropic body. 

(3.1) Homogeneous Deformations. 

This is a simple case where the strain tensor is constant throughout 
the volume of the body. In this case using (4.3) and the Lame relations 
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(1.6) and (1.7) we verify that the free energy per unit of volume F defined 
(3.1) becomes written as 

F = [E/2(l+v>] { u ik 2 + [v/(l-2v)] u (( 2 ) (3.4), 

o ik = [E/(l+v)] { u ik + [v/(l-2v)] u EE 8 ik } (3.5) 

and, conversely, 

u ik = [(l+v)oik-vo eE 8ik]/E (3.6). 

The total energy F total is obtained integrating (3.4) over volume V of the 
body, 

F to iai = Iv F dV (3.7). 

(4) Energy of Deformed Thin Plate. 

By a thin plate we mean that its thickness along z, normal to the 
plane (x,y), is small compared with its dimensions in the other two (x,y). 
The deformations themselves are supposed small, as before. In the present 
case the deformation is small if the displacements of points in the plate are 
small compared with its thickness. 

Let us suppose that the displacement vector u for points in neutral 
surface (middle of the plate) is given by 

u x — u y = 0 and u z = ^(x,y) (4.1). 

If there are only internal forces, that is, F e xtemai = 0 from (1.4) we have 
da ik /dx k = Oand the boundary condition n k o ik = 0, where n is unit vector 
outward normal to the film surface. [3] Since the plate is only slightly 
deformed we can take the normal vector n along the z-axis. Thus, in both 
surfaces of the plate o xz = o yz = o zz = 0. As the plate is thin and they are 
zero on each surface these tensor components must also be zero 
everywhere in the plate. We can therefore equate them to zero and use this 
condition to determine the components of the strain tensor. Thus, from the 
general formulae (3.5) we have 

Gzx —■ [E/(l+V)] U zx , O zy — [E/(l+V)] U z y 

and 


a zz = [E/(l+v)(l-2v)] { (l-v)u zz + v(u xx + U yy ) } (4.2). 
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Equating these expressions to zero, [3] we obtain cujdz = - du z /dx, 
duy/dz = - cu jcy, u zz = -v(u xx + u yy )/(l-v). In the first two of these equations 
the component u z can, with good accuracy, be replaced by £(x,y): 

dujdz = - dC/dx, ou y /oz = - dC/d y, whence 

u x = -z (dC/dx) and u y = -z(d'CJdy) (4.3). 

If the function £ = £(x,y) is known integrating (4.3) using the boundary 
conditions u x = u y = 0 for z = 0 we can determine u x (x,y) and u y (x,y). 
In the next step, knowing u x (x,y) and u y (x,y) we can calculate all 
components of the strain tensor u ik taking into account that 

u xx = -z(d 2 (Jdx 2 ) , u yy = -z(d 2 (Jdy 2 ) , u xy = - z{d 2 (Jdxdy), 

u xz =u yz =0 and u zz = z [v/(l-v)] {d 2 (Jdx~ + d 2 C/dy 2 ) (4.4). 

Now we can calculate the free energy per unit of volume of the plate 
using the general formula (3.4). A simple calculation gives 

F = z 2 [E/(l+v)] {(d\/dx 2 + d 2 (Jdy 2 ) 2 + [l/2(l-v)][(d 2 ^/3x3y) 2 - (d 2 (Jdx 2 ){ d\/dy 2 )]} (4.5) 

The total free energy of the plate F p i ate is obtained by integrating over the 
volume of the plate. The integration over z is from -t/2 to +t/2, where t is 
the plate thickness, and (x,y) over the surface of the plate. So, the total 
energy F p i ate = |F dV= jFdxdydz of the deformed plate is, with dA = dx dy, 

Fpiate = 

[Et 3 /24( 1 -v 2 )]JJ[(<9 2 ^/<9 x 2 + d 2 C/dy 2 ) 2 + 2(\-y){{d 2 C/dxdyf- (d 2 (Jdx 2 ){ d\/dy 2 )}]dA (4.5). 


(5) One Dimensional Thin Film Wrinkling. 

In Fig.2 a rigid film (f) is deposited on a soft elastic substrate(s). 
This rigid film attached to a soft elastic medium creates a large number of 
waves, wrinkles, on the system film & substrate. This wrinkling is a result 
of a compression of the film(see Fig.2). 
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Figura 2. Strain release of a film (f) connected to the soft substratum(s): in the upper 
side we see the film before the release. In lower side is seen after the release when the 
wrinkled state emerges with wavelength A and an amplitude C Q . 

(5.1) Elastic Energy of the Deformed Film. 

To estimate the elastic free energy F f of the wrinkled rigid film with 
thickness t we suppose that the film deformation £(r) = u(x,y), in the 
z-direction is given by 

£(r) = Co cos(kx) (5.1), 

where k = 2 tt/A, are pure cosine deformation in the x-direction and where 
« t. In this way (4.5) becomes written as 

F f =[Eft 3 /24(l-v f 2 )] J d 2 r |lapl(i;(r))l 2 (5.2), 

where E f and v f are, the Young modulus E f and Poisson ratio of the film, 
respectively. [3,6] The integral is done over the (x,y) plane of the film; the 
deformation is such that no stretching of the film is required. [3,6] The 
amplitude of the deformation is also assumed to be much smaller that the 
wavelength A, that is, A » C, Q . Inserting (5.1) in (5.2) and performing the 
integration over the (x,y) plane we get 

F t /A = [Eft 3 /48(1- v f 2 )] 2 k 4 (5.3), 

where A is film area. In Fig.2 is shown the film attached to a substrate (s) 
deformed by the wrinkle of the film (f). 

(5.2) Elastic Energy of the Deformed Substrate. 

The elastic energy F s of the deformed substrate in terms of the 
deformations u(x,y,z) is given by (3.7): [3,6] 
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F s = J v d 3 r [(.is (Ujk - 5 ik u u /3) 2 + K s u ££ 2 /2] 


(5.4), 


where the integration is done in the substrate volume V. Of interest is the 
case of low shear modulus p s = G s compared to bulk modulus K s = E s , that 
is, K s » p s or, equivalently a Poisson ratio v ~ 1/2 typical of elastomers. 

In this case using (3.5) we have a ik = K s u ££ 8 ik showing that only diagonal 
terms oi k contribute, that is, ou = 3K s u,j. As there is no external forces 
from (3.3) we have o xx = o yy = g zz = 0, resulting u ££ = 0, that is, 

div(u) = 0 (5.5). 

This equation must be solved with respect to the boundary conditions 
u(x, z = 0) = [0,0, ^ cos(kx)], i.e. the surface of the medium must match 
the wrinkle of the plane given by (5.1). Very far from the film the 
deformation must be zero: u(x, z — >oo ) = 0. Solving (5.5) submitted to the 
above boundary conditions we have 

u(x,y,z) = [ kzsin(kx), 0 , (l+kz)cos(kx)] ^ 0 exp(-kz) (5.6). 

Inserting (5.6) in (5.4) we obtain [6] 

Fs/A = G s k 02 (5.7). 

The wrinkle has emerged in the film due to the interaction between 
the film and substrate. We add the plate energy given by (5.3) with film 
energy given by (5.7) to obtain the total energy F tot of the system: 

F t0 , = F f /A + F s /A = [E,t 3 /48(1- v, 2 )] i^k 4 + G s k 02 (5.8). 

It is clear that there is an optimum wavelength X that minimizes the total 
F tot free energy. Farge wavelengths (low k) are not favorable due to the 
large deformations of the substrate, whereas the short X are too costly due 
to the bending of a rigid film. [6] From minimization of (5.8) with respect to 
k it is found that the optimum (1-direction wavelength) k ld is given by 

X,d = 2ntr\ id ,B (5.9), 

where 

ilid = E,/12G s (l-v f 2 ) (5.10). 
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Taking into account the Lame relation (1.7) G = E/2(l+v) we get 

il u = [(1 ■+ v s )/6 (1 - v, 2 )] (E, /E s ) (5.11) 

Thus, using (5.9) and (5.11) we verify that the 1-dimensional wrinkling 
wavelengths kid are given by 

k ld = 27it C sf (E f /E s ) 1/3 (5.12), 

where the parameter C sf = [(1+ v s )/6 (1- v f 2 )] 1/3 . 

Let us take, for instance, t = t* 1CT 9 m, k = k* 10 6 m, E f = e f 10 9 Pa 
and E s = c s 10 6 Pa = e s MPa. In this way, defining R = ( k*/27it*C s f) we 
get from (5.12) 

e s = 10' 2 e f / R 3 (5.13). 

A simple estimation of (5.13) can be done putting v s ~ v f ~ 0.3, 
t* ~ 2 , k* ~ 5 getting R 3 ~ 0.3 and, consequently, e s ~ 3.33 10 2 e f . In this 
way, if E f ~ 80 GPa, that is, ef = 80 we see that e s ~ 2.6, which implies that 
E s ~ 2.4 MPa. This value is compatible with measured Young polyurethane 
elastic modulus. [7] 

(6) Isotropic Thin Film Wrinkling. 

Supposing that instead of a 1-dimensional wrinkling, analyzed 
above, there is an isotropic wrinkling, the parameter rj ld would be replaced 

by r| iso : [6] 

r|iso= {(3-4v s )/(l- v s )} q ld (6.1), 

In this way instead of (5.9) we have now 

k iso = 27Tt Csf(iso) [E f /E s ] 1/3 (6.2), 

where C st (iso) = [(3 - 4v s ) (1+ v s )/6(l- v s )(l- v f 2 )] 1/3 . 


(7) Experimental and Theoretical Results. 

In this section our theoretical prediction are compared with the 
experimental results of Fernanda et al. [8] shown in Fig.3.In this experiment 
the wrinkling wavelengths have been measured when DLC films are 
deposited on PDMS substrates previously exposed to oxygen plasma. 
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DLC 10 nm 


DLC 20 nm 


Figure 3. Wrinkling wavelengths measured in reference [8]. 

The measured wavelengths X for t = 19.5 nm with E f = 75.8 GPa are, 
Mt) = 2.200, 3.460, 5.870 and 11.680. (nm) (7.1) 
and for t = 10.8 nm with E f = 192 GPa are 

A,(t) = 1.660, 2.600, 3.520 and 4.560. (nm) (7.2). 
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(7.1) One-dimensional Wrinkling. 

For t = 19.5 nm we verify, using (5.12) and (5.13), that E s measured 
in MPa would be given, respectively, by 4.4, 1.1, 0.2 and 0.03. Showing a 
fair agreement between theory and experiment. For t = 10.8 nm we verify, 
using (5.12 and (5.13), that E s measured in MPa would be given, 
respectively, by 4.4, 1.1, 0.4 and 0.2. Showing a fair agreement between 
theory and experiment. 

(7.2) Isotropic Wrinkling. 

Putting v s ~ Vf ~ 0.3 we verify that C S f(iso) ~1.6 C S f(ld) implying 
that k iso ~ 1.6 kid .This shows that the wavelengths of isotropic wrinkles 
would be about 2 times larger than the 1-dimensional ones. 
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